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Abstract 

We investigate a D dimensional generalization of the Schroedinger- 
Newton equations, which purport to describe quantum state reduction as 
resulting from gravitational effects. For a single particle, the system is a 
combination of the Schroedinger and Poisson equations modified so that 
the probability density of the wavefunction is the source of the potential 
in the Schroedinger equation. For spherically symmetric wavefunctions, a 
discrete set of energy eigenvalue solutions emerges for dimensions D < 6, 
accumulating at D = 6. Invoking Heisenberg's uncertainty principle to 
assign timescales of collapse correspoding to each energy eigenvalue, we 
find that these timescales may vary by many orders of magnitude depend- 
ing on dimension. For example, the time taken for the wavefunction of a 
free neutron in a spherically symmetric state to collapse is many orders of 
magnitude longer than the age of the universe, whereas for one confined 
to a box of picometer-sized cross-sectional dimensions the collapse time is 
about two weeks. 



1 email: rgmelko@gandalf.uwaterloo.ca 
2 email: mann@avatar.uwaterloo.ca 



1 Introduction 



In quantum mechanics, objects are described by wavefunctions. These take the 
form of complex superpositions of various evolutionary alternatives, or states. 
Although successful in describing many aspects of the quantum world, this pic- 
ture often leads to troubling interpretations when extrapolated to the macro- 
scopic level. One issue that has suffered long debate is the fact that one never 
observes a superposition of states. Rather, one only observes a system's basic 
or stationary states. We are therefore forced to provide a mechanism by which 
quantum wavefunctions reduce to their stationary states. This process is called 
wavefunction collapse or state reduction. Motivated by the basic conflicts which 
exist between general relativity and quantum mechanics, a number of authors 
have proposed the idea that wavefunction collapse is an objective phenomenon 
which arises due to gravitational effects jij. For example Penrose has sug- 
gested a scheme in which a superposition of two stationary quantum states 
should be fundamentally unstable if there exists a significant mass displacement 
between them. In this case there should be some characteristic timescale Tq for 
decay into the basic states. Although a detailed estimate of Tq would require a 
full theory of quantum gravity, under this hypothesis it is reasonable to expect 
that for non-relativistic systems 

where AEq is the gravitational self-energy of the difference between the mass 
distributions of the two states. 

The explicit nature of the basic states in this consideration is somewhat 
unclear. We cannot simply regard the position of a lump of mass as a basic 
state, because then we would be forced to regard any general state of a particle 
as a superposition. As a possible solution to this problem, Penrose proposes that 
these (non-relativistic) basic states are solutions of the Schroedinger equation 

V 2 * + U^ = E^ (2) 

2m 

where the additional term represents a coupling to a certain gravitational poten- 
tial U. This potential is determined (via the Poisson equation) by the expecta- 
tion value of the mass distribution in the state determined by the wavefunction. 
For single particle systems, the matter density is determined by the probability 
density from the wavefunction, and so 

V 2 U = ATrGm 2 |*| 2 (3) 

where G is Newton's gravitational constant, and m is the mass of the single 
particle. Equations (^,||) are dubbed the Schroedinger-Newton (SN) equations 
[pi. A preliminary investigation of the properties of the solutions to the SN 
equations was recently carried out by Moroz et. al. Q|. Under the assumptions 
of spherical symmetry in 3 dimensions, and by demanding only that U and $ be 
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everywhere smooth, they discovered a discrete family of bound state solutions, 
labelled by an integer n > 0. Each solution is a normalizable wavefunction, 
and the nth solution has n zeros. The energy eigenvalues associated with each 
of these solutions are negative, and monotonically converge to zero for large 
n. These results can be justified analytically 0. The energy eigenvalues are 
the differences between a given bound state and a continuum 'superposition' 
state, and so provide via (Q) an estimate of the timescale of self-collapse of a 
single particle of mass m. The energy eigenvalues scale like m 5 , and so particles 
of small mass have extremely long self-collapse times - for a nucleon mass the 
estimate is 10 53 s B. A recent related study by Soni is commensurate with these 
results ||. 

Relaxing the assumption of spherical symmetry is in general a difficult task 
due to the non-linearity of the SN equations. However there are two situations 
in which this is fairly straightforward: cylindrical symmetry with no angu- 
lar momentum and planar symmetry. Rewriting the SN equations for these 
cases effectively reduces them to 2 and 1 dimensional situations respectively. 
These cases, along with the spherically symmetric case, can be simultaneously 
recovered by rewriting the spherically symmetric SN equations in D dimen- 
sions. Motivated by the above, we consider in this paper an analysis of the 
-D-dimensional spherically symmetric SN equations, for D > 1. Although the 
higher-dimensional cases are of less direct physical interest that the D = 2, 3 
cases, such a study affords us some insight into the dimensional behaviour of 
the SN system. This behaviour may be of more than pure pedagogical interest 
since many candidate approaches to quantum gravity are typically cast in higher 
dimensions (superstring theory being the obvious example). 



2 The D-dimensional SN equations 

Any solution to the SN equations (|],||) must be normalizable (i.e. square- 
integrable). Integrating the probability density over all space yields 

/>oo 

k 2 = d D x |*| 2 (4) 
Jo 

where A: is a dimensionless number, and so the wavefunction must be rescaled 
to ensure there is unit probability of finding the particle somewhere in space. 
Writing ^ = ktp, the SN equations then become 

Eil) (5) 
AitGk 2 m 2 |i/>| 2 (6) 
1 (7) 



2m ^ 



V 2 U 



d D x ivr 
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and we see that the normalization factor enters the system due to its non- 
linearity. Redefining variables in dqj9) via H] 



where 



ip = aS E-U = (3V (8) 



yields 



V8irGk 2 m 3 k 2m 

V 2 5 = -SV (fO) 
V 2 V = -S 2 (ff) 



a 2 



d D xS 2 = 1 (12) 



where we can assume that ip is real without loss of generality. The parameters 
a and (3 have units of (length) 2 D ^ 2 and (length) 2 x energy respectively. The 
system (|l0 11) is invariant under the rescaling transformation 

(S, V, x) — ► (A 2 ^ X 2 V, X^x) (13) 

independent of the dimension D, where A has units of inverse length. Using this 
transformation we can rewrite the system in terms of fully dimensionless func- 
tions (S, V) of dimensionless variables. For the spherically symmetric case the 
D-dimensional Laplacian operator is V 2 (p = r l ~ D (r " 1 ^')' and so ( |To| , |l~T] , |T2| ) 
become 

(r^S')' = -r D ~ l SV (14) 
(r D -V)' = -r^S 2 (15) 

poo 

a 2 X^ D \ d D xS 2 = 1 (16) 







where the prime denotes differentiation with respect to r. We must require that 
the integral in (|l^) be finite in order for the SN system to be physically meaning- 
ful. We therefore seek solutions to the SN system that are finite for all values 
of r > 0; continuity requirements imply that these solutions are everywhere 
smooth. Note that a^X 40 is a dimensionless quantity. 

However, our equation (|l6| ) clearly has a problem for D = 4, where we lose 
our ability to rescale using A. Because we still require the wavefunction to be 
normalized to unity, we are forced to introduce a constant & > in the D = 4 
case so that S(r) — > V&S(r). Then our equations (Q), (|l5|) , and ( |l6| ) become 



{r 3 S')' = -r 3 SV (17) 

(r 3 V)' = -&r 3 S 2 (18) 

O 

a 2 £ / d 4 xS 2 = 1 (19) 
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We will not always include this constant & explicitly in the following discussion. 
It will simply be understood to occur in the rescaling when D = 4. Smoothness 
implies that S'(0) and V'(0) are finite. With this information we can rewrite 
0n>|) as 

S(r) = So + j^J\(l-(^) D ~ 2 ^S(x)V(x)dx (20) 

V(r) = V + j^J\(l-^) D ^S 2 (x)d x (21) 

It is straightforward to show that Picard's theorem Q is satisfied by this system 
of equations, and so given So and Vq a unique solution to ( |l^JT5| ) exists within 
a range [0, R (So, Vq) ). The D < 3 versions of the integral equations ( p(|^l| ) 
require some care. Integrating (jl4||l5|) explicitly for these cases yields 



S(r) = S Q + (x-r)S(x)V(x)dx (22) 
Jo 

V(r) = V + I (x-r)S 2 (x)dx (23) 



for D — 1 and 



S(r) = So + j xln(~^)S(x)V(x)dx (24) 
V(r) = V + f xln(-) S 2 (x)dx (25) 



for D = 2 . Both of these sets of equations may be obtained formally from 
( pO 21) by inserting these values of D, the latter case being understood as the 
limit D -> 2. 

Note that for all of these cases the function V(r) is montonically decreasing 
since 

f r fr\ D-i 

V'(r) = -J {=■) S 2 (x)dx (26) 

a result valid for all D > 1. Hence if Vo < 0, then limr^oo V(r) = — oo and 
so S(r) will also diverge for large r. Consequently only Vq > is of physical 
interest. For Vq > the function V(r) > initially. By rewriting S(r) — 
r~%- s(r) equation ( |l4|) can be rewritten in the form 



r 2 s" + rs' + 



r 2 V(r) - 



s = (27) 



which for V(r) = Vo is Bessel's equation. Hence near the origin we expect s(r) 
to have oscillatory behaviour. However when V(r) becomes negative eq. ( p7| ) is 
similar to the modified Bessel equation and the behaviour of s(r) will be a linear 
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combination of exponentially amplified and damped functions. Normalization 
requirements imply that only the exponentially damped solutions are allowed. 
In this case S(r) decays like an exponential times an inverse power of r and 
so the integral in ( pl| ) will be finite, yielding a finite V(r) at large r. This 
behaviour was already noted for the D = 3 case [Q; we see here that it is valid 
for arbitrary D > 0. 



3 Numerical Study of the D-dimensional SN Equa- 
tions 

The form ( |l4| , |l5| ) of the SN equations is not well-suited for numerical study. A 
more convenient form is 

(rS)" = -rSV + (3- D) S' (28) 
(rV)" = -rS 2 + (3- D)V (29) 



We are interested in finding solutions to the system (|2§|,|29j) for that are smooth 
and finite for all r. These are the bound state solutions referred to earlier. We 
will require our initial conditions, 5(0) = So and V(0) = Vq to be greater than 
zero, a constraint which avoids both the trivial solutions S — V = and the 
non-normalizable solutions V — ±S = —2(D — 4)r~ 2 , and is consistent with 
the rescaling freedom Jig). It is straightforward to obtain the power-series 
solutions to (plJ29) near the origin: 



2\ 



af \ c S oVo 2 S (S 2 + V 2 ) 4 S V ((5D + 4)S 2 + DV Q <j 6 

b{r > ~ bQ ~^D 7 ' + 8D(D + 2) r " 48£>2(D + 2)(£> + 4) T W 

[) ° 2D AD(D + 2) 24D 2 (D + 2) (D + 4) [ ' 

where the functions S and V and their first derivatives are required to be finite 
at the origin. Using the rescaling freedom to set Vq — 1 (which is equivalent to 
setting A 2 = Vq), we see that for any dimension D the solutions near the origin 
depend only on a single free parameter. 

Integrating the system (|2^j2^) using a Fehlberg fourth-fifth order Runge- 
Kutta method in MAPLE we find that in each dimension D an infinite set 
of discrete bound states appears as expected. These can be identified by the 
number of local extrema. An infinite amount of precision is required in the 
choice of So such that the solutions do not diverge as r increases. This value of 
So marks the transition between solutions in which S(r) diverges to +oo and 
solutions where S(r) diverges to — oo. As we increase the accuracy of So for a 
specific bound state, the value of r at which S(r) blows up increases. We will see 
that the values for So and the distance between bound states vary significantly 
with dimension. 

The general method for finding solutions to the D-dimensional system is 
the same as that in the D = 3 case [Q. Choose a value of So = Sq™ +1 ' + 
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for which S(r) has n zeros and diverges to (— l) n oo at finite r. Then select 
a slightly smaller value of Sq = Sq" +1 ^ which has n + 1 zeros and diverges 
to (— l)™ +1 oo for some finite r. Between these two values of Sq is a value 
Sq +1 '~ < 5q™ +1 ^ < S'q™ +1 ' )+ for which S(r) has n zeros and is smooth and finite 
for all values of r and is square-integrable. The potential V(r) will also be 
smooth and finite, with finite energy eigenvalue E n+ \. By successively narrowing 
this interval the bound-state value of Sq can be achieved to any desired accuracy. 
The bound state wavefunction with n + 1 zeros may be obtained by choosing 
another value S u ™ +2 '' + < So" 1 ' - and repeating the procedure with n — > n+ 1. 
We find that the values Sq +1 ' are rapidly decreasing functions of the dimension 
D, and that the bound state solutions appear to accumulate at D = 6. We are 
limited by numerical accuracy in carrying out our investigations for dimensions 
D > 6 due to this effect. A summary of this data is presented in Table 1. 



Si 



S' 2 
a 



St 



S' 1 



1.5583977884 
1.2134344293 
1.0886370794 
1.0327684253 
1.0081005592 
1.0000000000 
Initial conditions Sr 



0.379904201 

0.6482524055 

0.8264742841 

0.9309542414 

0.9825730584 

1.0000000000 

n+l) 



0.2128374651 0.1475990005 

0.4937184140 0.41447908 

0.7442133785 0.70014479 

0.9053504436 0.89422924 
0.9795819080 > 0.978 
1.0000000000 1 

Table 1. Initial conditions Sq 1 ^ 1 ' for the first (n — 0) to fourth (n = 3) bound state 
wavefunctions, for one to six spatial dimensions. The accumulation of the higher 
bound states for large dimensions results in the decreasing accuracy of the lower 
right hand entries of the table. 

Figures 1,2 and 3 illustrate the solutions for S for D = 1, for the first, second 
and third bound state transitions. Figures 4 to 8 illustrate the first bound state 
transitions for dimensions D — 2 to D — 6. Figure 9 illustrates the accumulation 
effect at D = 6. 



We turn now to the problem of determining the energy eigenvalues for the 
bound states. From (fL5|) and (^l|) we find that the large-r expansion of V(r) is 



V(r) = A D + 



where 



Ad 



Bd = 



V 



Bp 1 
D - 2 r D " 2 



2-D 



x D - 1 S 2 (x)dx 



xS 2 {x)dx 



(32) 

(33) 
(34) 
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Figure 1: First bound state transition for the ID SN equations 
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Figure 2: Second bound state transition for the ID SN equations 
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Figure 3: Third bound state transition for the ID SN equations 
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Figure 4: First bound state transition for the 2D SN equations 
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Figure 5: First bound state transition for the 3D SN equations 



m 



-1- 



-2 



So = 1.03276845 



So = 1.03276840 



10 20 30 40 50 



Figure 6: First bound state transition for the 4D SN equations 
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Figure 7: First bound state transition for the 5D SN equations 
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Figure 8: Integration of the 6D SN equations. Bound state characterictics at 
this and higher dimensions are not apparent. 
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Figure 9: S(0) versus dimension. Bound states appear to accumulate at D = 6. 



These expansions are formally valid for all D > 0; the explicit expansions for 
D < 3 are 

f Ax -B 1 r + --- D = l 

V ^ = \ M-B,ln(s-) + ... D = 2 ^ 



and 



V + fi°xS 2 (x)dx D = l 



Ad - { v„ f°°~ir, fx Vc2a/m~ D = 2 ^ 



Bn I £ S2 ^ dx D = 1 (37) 

where ro is an arbitrary length scale. The energy eigenvalue is given by 

E D =PV{r D )=PA D (38) 

The quantity \d is the point at which the potential U(r) vanishes. For D > 3, 
Xd = oo. However for D — 1,2 the potential diverges there and the situation 
is more delicate. For D = 1 a Newtonian gravitational potential diverges 
linearly with r = |x| . Here we extract the r-independent term from ( pq) to 



obtain Ed=\] the result is still given by the right-hand side of (42) with D = 1 



Effectively we have set xd — in (|35|) even though this equation is a large r 
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expansion. For D — 2 the Newtonian potential diverges at both large and 
small r. The normalization point xd = Xq is therefore arbitrary, and we set 
Ed-i = (3X d=2 Ad=2 so that (E^) remains valid in this case as well. This is 
tantamount to setting U(xq) = J r ^°xln (—^ S 2 (x)dx . A natural normalization 
point is to choose Xq to be the point at which the potential V(xq) = 0. Note 



that under the rescaling transformation ([13]) (A, B) — ► ( X 2 A, A 4 D B ) , where 



A D = ^ + ^—J^j xS 2 (x)dx (39) 



(i 



oo 

B D = I x D - 1 S 2 (x)dx (40) 



and so the combination A D D /B 2 D = A D D / B D is invariant. The parameter A 
(or alternatively Vo) is completely arbitrary, serving the sole function of setting 
the length scales of the problem in units of aj^ 4 D ^ ■ Solving (|l^) for A yields 



r(f) 



2ir 2 cx 2 Bd 
With this normalization, we find that 



(41) 



/ r (—) \ i ~ D a 

E D = PA D = P\ 2 D A D = 13 — W- ?s- (42) 

\l-K-a 2 ) (Bd)^ 

and so the D-dimensional energy is given in units of /3/a 4 /( 4 --°) . The quan- 
tities Ad and Bd are straightforwardly solved numerically from ([53],[34]). The 
preceding expression can be rewritten as 

E D = -m pl c 2 [ y2 _{ ^- (43) 

2 \m p ij y n — J (B D )~ 

where m p i = (h D ~ 2 j 'G j)C D_4 ) 1 " D is the Planck mass in D > 1 dimensions 
(in one spatial dimension the quantity c 3 /HGd is unitless, and can be absorbed 
into the normalization constant k). 

For D < 3 the energy eigenvalue is an increasing function of the particle 
mass, whereas for D > 5 it is a decreasing function. The most rapid increase is 
Ed ~ m 5 for D = 3, with only linear and quadratic dependence in D = 1 and 
2 respectively. The time-scale for collapse of the state vector is therefore most 
rapid in a world with three spatial dimensions for any bodies whose mass is ap- 
preciably larger than the Planck mass. In a world of more than four dimensions 
the collapse is fastest for the lightest-mass particles, leading to behaviour which 
is at complete odds with that expected in the macroscopic world. 
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If we interpret the lower-dimensional cases to be situations in three spatial 
dimensions with planar or cylindrical symmetry, a somewhat different picture 
emerges. Inserting the constants into (^3j), the dependence of the energy eigen- 
values for the planar, cylindrical and spherically symmetric cases are respectively 

Ep = -^ mpl c — J -) — (44) 



2 



pU \ J (Si) 
LA A 



5 

2 



E s = 2m pl c' — ^ (46) 

\m p ij (B 3 f 

where £ p i is the Planck length and £ is the dimension of the box (or cylinder) in 
which the particle is confined. From the perspective of state-vector reduction, 
a particle prepared in a plane-wave or cylindrical state will have a radically 
different time-scale of collapse than the same particle prepared in a spherically 
symmetric state. A neutron in a spherically symmetric state will have a collapse 
time of At ~ 10 53 s whereas a neutron confined to a square box 10m in cross- 
sectional size will have a collapse time of only At ~ 10 22 s. Both of these times 
are much longer than the age of the universe. However by confining a neutron 
to a rectangular pipe whose width is on the order of 10 picometers, the collapse 
time becomes on the order of 10 6 s, or about 12 days. Such an experiment 
would push the limits of current technologly. 

Note that for D — 4 the probability is scale- invariant, and it is not possible 
to solve ( fl6| ) for A. The energy scale in this dimension is independent of the 
wavefunction normalization. If we work through the problem using the & 
rescaled equations ( |l8| ) and (|l9|), we can recover an expression for the energy 
eigenvalues in this dimension using (^) . When rescaling equations ( |32] ) we can 
choose to absorb ^ in the rescaling of the expansion coefficients: (A, B) — > 
(AX 2 A, RB) where A and B are equivalent to (||) and (j40j) . Then solving 
for & we get 



R = T ^ (47) 
We use this normalization to solve for the energy eigenvalue when D = 4 

E 4 = (3A 4 = f3&\lA 4 =\ 2 4 (3 (^Pl) IT ( 48 ) 

\2-K z a z J B4 

Clearly we are left with an unsolved parameter, A4, which is at this point com- 
pletely arbitrary. This is a consequence of the introduction of the rescaling 
constant & in the D = 4 case. 



13 



A graph of the invariant A^ 4 ~ D ^ / B 2 is instructive because it reflects the be- 
haviour of the energy eigenvalue solutions with bound state that is independent 
of the value of the D-dimensional Newton constant. Figure 10 shows that the 
energy eigenvalues are positive and increase with bound state for D < 3, consis- 
tent with what we expect for the form of the classical Newtonian gravitational 
potential in these dimensionalities. Figure 11 shows that the energy eigenvalues 
are negative for dimensions greater than 3. In this graph, the values of A4 and 
r are set to unity for D = 4. 

4 Conclusions 

We have presented the preliminary results of our numerical and analytical in- 
vestigation of the Schrocdinger-Newton equations in D dimensions under the 
assumption of spherical symmetry. Consistent with previous analyses in 3 di- 
mensions, we have demonstrated numerically the existence of a discrete set of 
"bound-state" solutions for the equations, which are associated with different 
energy eigenvalues. The bound states appear to accumulate for dimensions 
greater than or equal to six. The bound state energies for D < 3 are positive 
and increase with increasing bound state number. For D > 3 the energies are 
negative and converge to zero with increasing eigenvalue number. This is consis- 
tent with Moroz Q , who found that for D — 3 the energy eigenvalues converge 
to zero as the order of the bound state increases. For higher dimensions, the 
variation of the eigenvalues with bound state becomes less pronounced, as illus- 
trated in Figure 11. For D > 5 we find that the eigenvalues accumulate, and 
we are unable to determine any energy differences between bound states to 10 
significant digits. 

Moroz ventured further in the three dimensional study to assign numerical 
estimations for timescales of collapse, using the lowest bound state solutions to 
determine energy eigenvalues. However the analogous calculation using general 
D dimensions is not strictly correct due to possible dimensional variation in 
Newton's constant Gd which appears in (g). We will not attempt to predict 
the variation of Gd with dimension here, although we note that Kaluza-Klein 
theory suggests that Gd is proportional to the three dimensional G3 divided by 
the square root of the volume element of the extra space dimensions. We leave 
any such numerical estimates for future work. 

However if we interpret the lower-dimensional cases to be systems in three 
spatial dimensions with either planar or cylindrical symmetry, a different picture 
emerges. The collapse time scale depends not only on the mass of the particle 
but also on the size of the box in which it is confined. If the gravitational 
influence of the box can be neglected and if the gravitational field of the neutron 
in a plane wave state has planar symmetry, then we find that the collapse 
timescales differ enormously between the two situations. A free neutron under 
its own gravitational self-influence will experience a collapse of its wavefunction 
on a timescale many orders of magnitude larger than the age of the universe, 
whereas one confined to a rectangular pipe whose cross-sectional width is about 
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Figure 10: Behaviour of the invariant A 1 - 4 D ' / B 2 for the first three bound states 
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10 picomctcrs in size will decay in about two weeks. 

An obvious extension of our work involves introducing external field into our 
version of the SN equations, in the form of a point source perturbation. This 
can be added to the Schroedinger equation to give the perturbed system 

fi 2 M 

V 2 * + C7* * = (49) 

2m r 

V 2 U = 47rGm 2 |*| 2 (50) 

Here M is the mass which generates the perturbing field. This system is quite 
interesting in the respect that the form of the bound state solutions is drasti- 
cally altered from that of the "free field" equations considered in this paper. A 
preliminary study of this system indicates that each bound state solution seems 
to have three regions. For small M, the solutions asymptotically approach those 
of the free field or non-perturbed case, as expected. For moderate values of M, 
the equations are goverened by both the free field and point source field, causing 
the values of So, for which the solutions do not diverge, to alter considerably. Fi- 
nally, for very large M, the solutions are goverened almost entirely by the point 
source perturbation, and become similar to the hydrogen atom solutions of the 
Schroedinger equation. The effect of incrcsing dimension causes the "transition" 
values of M (which separate the different regions of the solutions) to decrease. 
A full study of the point source perturbation is currently in progress. 
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